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Abstract 

We propose a "multi-stream" inflation model, which is a double field model with spontaneous 
breaking and restoration of an approximate symmetry. We calculate the density perturbation and 
non-Gaussianity in this model. We find that this model can have large, scale dependent, and 
probably oscillating non-Gaussianity. We also note that our model can produce features in the 
CMB power spectrum and hemispherical power asymmetry. 
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I. INTRODUCTION 



Inflation has become one of tlie most important windows to connect cosmological ob- 
servations to fundamental physics, such as string theory. Recent studies on string theory 
show evidence that inflation may happen in the string theory landscape [l| with extremely 
complicated inflaton potential 2|. 

One lesson from the string landscape is that the inflationary dynamics may be not as 
simple as one expected before, and more possibilities during inflation should be taken into 
consideration. Inflation in the string landscape opens up a number of possibilities to test 
string theory by the stringy imprints on the CMB temperature fluctuations. 

On the other hand, in the recent analysis of CMB anisotropy, there are several hints 
for new physics. If future experiments confirm some of these possibilities, they will pose 
challenges for inflationary cosmology and opportunities for string theory. These challenges 
and opportunities include: 



Non-Gaussianity. The WMAP5 bound for local non-Gaussianity is 



9 < /^i < 111 (95%CL) . 



(1) 



The result is still consistent with f^L — 0. However, the central value of /nl is large 
and positive. If integration of more years of WMAP or the Planck mission confirm this 
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17l |. As there is a 



consistency relation 13|, |l8| for single field inflation, a large and local non-Gaussianity 
shows evidence that inflation should involve more than one field. 



Features in the CMB. As shown in 19|, there may be features in the CMB scalar power 
spectrum. A feature in the CMB spectrum implies a departure from the standard single 
field slow roll inflation at a certain scale. 

n 

Hemispherical power asymmetry. As shown in |[20|], the temperature- fluctuation am- 
plitude is larger by roughly 10% in one hemisphere of the CMB map than in the other. 
Again, this asymmetry can not be explained by the simplest inflation model. It would 
also be interesting to see whether this difference of power also happens on smaller 
scales in the sky. 



In this paper, we propose a "multi-stream" inflation model, which is a double field infia- 
tion model with spontaneous breaking and restoration of an approximate symmetry. We find 
that with appropriate choice of parameters, our model is consistent with current observa- 
tions, and can produce large non-Gaussianity, features in the CMB and hemispherical power 
asymmetry. In Section 2, we describe the picture of multi-stream infiation and calculate the 
density perturbations. In Section 3, we calculate the non-Gaussianity of multi-stream infia- 
tion. In Section 4, we give a summary of the signatures for multi-stream infiation in different 
parameter regions. We conclude in Section 5. 

II. MULTI-STREAM INFLATION 

In this section, we investigate the background evolution and the scalar perturbation 
spectrum of multi-stream inflation. We consider a period of inflation involving two flelds. 

We consider the potential V{ip, x) as illustrated in FiglH In the beginning, the ip direction 
is the inflationary direction. As ip evolves to (pi, the x direction becomes a tachyonic 
direction. At this point, the inflationary trajectory spontaneously breaks into two paths, 
namely A and B in Fig{TJ The whole picture looks like a stream which splits and flows along 
both sides of a hill. So we call this scenario "multi-stream" inflation. 

As inflation continues, depending on the potential and the reheating time, the trajectories 
A and B may either recombine into a single trajectory (as shown in FigH]) or not. As we 
will see in the calculation of scalar perturbations, if A and B do not recombine and reheat 
separately, one typically has large anisotropy at the symmetry breaking scale, and requires 
more flne tuning to flt the CMB power spectrum. So in this note, we mainly discuss the 
former possibility. 

Now let us investigate the perturbations. The perturbations of multi-stream inflation 
should be the same as in usual inflation models except in the regime between ipi and <^2- 
Moreover, if the potential energy of path A and path B are identical, multi-stream inflation 
also does not differ much from usual inflation models except for some amount of isocurvature 
perturbations, because in this case, to choose which way is not important. 

To have more interesting physics, we assume the paths A and B have a slightly different 
potential energy 6V = V{(p,xa) — V{'^,Xb), with 6V non-vanishing only if ipi < Lp < ip^- 
After this setup, multi-stream inflation will have several interesting differences from the 
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FIG. 1: In this figure, we illustrate the picture of multi-stream inflation. As the inflaton rolls 
down the potential, it experiences a spontaneously breaking of approximate symmetry at /ci, with 
if = ifi. After that, the inflation trajectory branches off into two paths denoted by A and B. From 
the viewpoint of the CMB map, some parts of the sky with size of order l/ki come from path 
A, and other parts of the sky come from B. At A;2, with ip = (p2, the approximate symmetry is 
restored and the inflationary direction becomes unique again. 

usual inflation model. 

Firstly, the perturbations along trajectories A and B have different power spectra. It is 
because the tilt of the potential is typically different when 5V 7^ 0. To be more explicit, the 
power spectra along the two trajectories are calculated using the local Hubble parameter 
and the local slow roll parameters as 



^ 87r2M2e 



(2) 

B 



The difference in the power spectra satisfies 

{P^ - Pf)/P^ ~ 6V/V . (3) 

In position space, trajectories A and B correspond to different parts of the sky. So 
a difference between the power spectra P^ and P^ implies that on cosmological scales 
corresponding to (fi and (f2, the perturbations in some regions of the sky have more power, 
and in some other regions of the sky have less power. 
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For example, if ki corresponds to the scale of the whole observable universe, then trajec- 
tories A and B both span about half of the sky. In this case, multi-stream inflation provides 
a possible solution to the hemispherical power asymmetry problem 20 1. 

If the hemispherical power asymmetry is really there and with a difference of order 10% 
in magnitude, we need 6V/V ~ 10~^, and ki to be the scale that comparable to the present 
Hubble scale to explain it. As shown in the following sections, in this case, the corresponding 
non-Gaussianity is typically /nl ~ O(IOO). 

In multi-stream inflation, ki could also correspond to much smaller length scales compared 
with the present Hubble scale, depending on the explicit form of the inflationary potential. 
In this case, the multi-stream effect generates about {ki/k^Y regions on the CMB, where fco 
is the comoving wave number whose wavelength is the whole observable universe. In these 
regions, the primordial power spectrum either takes the form or . 

Secondly, there can be a difference of e-folding numbers for paths A and B. This difference 
results in a time delay effect for the reheating surface, thus producing extra perturbation. 

To be explicit, let 5N = Na — Nb, where Na and A''^ are the background e-folding 
numbers evaluated along trajectories A and B respectively. We note that 5N depends on 
the detailed form of the potential 

where a prime denotes a derivative with respect to ip. We denote the comoving scale that 
exits the Hubble radius when (y9 = by fci, and = v^2 by k2 as shown in Fig{T] 5N provides 
an extra source of curvature perturbation on the scale ki. It is because the existence of 5N, 
thus the time delay in reheating divides the CMB sky into patches with comoving length 
scale 1/ki. This anisotropy corresponds to a perturbation 5C|^. The amount of perturbation 



can be calculated as 



2l| 

H 



Ck, =AN= -6^ + SN = Co + 5C , (5) 

where AA^ is the total e-folding number difference between path A and path B. The dif- 
ference AA^ has two sources, namely the quantum fluctuation of the inflation 6ip, and the 
intrinsic difference in the e-folding number SN between the two paths. 
The total power spectrum takes the form 

n.,-Co^+«)' • (6) 
5 



We neglected the cross term CoS(^^, because if one does not consider the conversion from 
isocurvature perturbation to curvature perturbation, then and 6(^^ are uncorrelated. We 
will return to this issue and give a more careful analysis of Eq. in the next section. 

The existence of 6N is a mixed blessing. The good news is that 6N provides a nice 
explanation for the features in the CMB power spectrum, because the jump of ( takes place 
only in a very sharp range of k. The bad news is that when 6V is not very small, to make 
this jump not too large, one needs some amount of tuning to be consistent with experiments. 
We shall estimate the amount of tuning in Section 4. 

III. NON-GAUSSIANITY OF MULTI-STREAM INFLATION 

In this section, we calculate the non-Gaussianity of multi-stream inflation. We shall 
neglect the non-Gaussianity produced by mechanisms other than the multi-stream effect. 

To perform the calculation, we further assume that the mass of x shortly before the 
Hubble exit of ki is of order H. This assumption guarantees that is not too light to 
randomly walk more than one e-fold away from the origin, and also guarantees that m^^ is 
not too heavy to have an effect. 

As the inflation trajectory turns at (fi, a certain amount of isocurvature perturbation 
6xki is projected to curvature perturbation. Note that the sign of 6xki decides whether 
inflation will be along trajectory A or B. The observed curvature perturbation at ki also 
includes a contribution from 6N, we have 



where c denotes the fraction of isocurvation which is projected onto the curvature pertur- 
bation during the change of the trajectory, 6(^_^ and 6(^_^ denote the curvature perturbation 
from projection of isocurvature perturbation and from 6N respectively. 

There are two kinds of non-Gaussianities in multi-stream inflation. Firstly, the one-point 
probability distribution function for (ki is non-Gaussian. This is because the extra density 
fluctuation S(^^, coming from the e-folding number difference between trajectories A and B, 
peaks at two particular values: 



Cfci = Co + ^Ck 



5Ck^ = c—5xk, +5N = 5Cl, + 5C 



■N 
ki 1 



(7) 



6N 



) 



(8) 



6 



Here we use P{x) to denote the probability distribution function. This should not be con- 
fused with the power spectrum P^. The contribution from the entropy perturbation 6(^_^ is 
still Gaussian. 

Note that here we assumed for simplicity that the probability to pick up trajectory A and 
B are equal. This seems natural in our discussion because we are talking about a spontaneous 
breaking of an approximate symmetry. One can generalize the discussion straightforwardly 
to the case with unequal probability between trajectories. With unequal probability, we can 



also provide a theory for the CMB cold spot in the multi-stream inflation framework j22|. 

Note that the one point distribution function is not very easy to probe, especially for 
one particular value ki. So we now proceed to discuss the other kind of non-Gaussianity 
produced by multi-stream inflation. The idea is that the amount of perturbation at comoving 
wave number ki and the amount of perturbation at comoving wave number ki < k < k2 are 
correlated. The correlation comes from the fact that the following three effects are 100% 
correlated with each other: 1. An isocurvature perturbation deciding whether to go along 
A or B. 2. The difference between and . 3. The difference of the e-folding number 
6N between A and B. 

The three point correlation function can be written as 

iCkACk) = (Coaa) + {sc^Mk) + • (9) 

As we have neglected non-Gaussianities from sources other than the multi-stream effect, we 
have 

{(oCkCk) = . (10) 

The two terms from entropy modes and SN can be written as an integration over the product 
of the three points and the probability distribution function, 

/oo 
dsc',fdap{sc',f, a)<5craa • (n) 
■00 

with the combined probability distribution function of SC^^^ and Cfc, 



>2 ^2 

e ^ e ^ 



(12) 
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given in Eq. with (erf) = ((CfJ ), <^a= (CDa and a% = {CI)b- is the Heaviside 
step function. To avoid 5 functions, we use the normahzation that the fluctuation is averaged 
over one Hubble volume. (A familiar example is that in this normalization, the fluctuation 



of the inflaton takes the form 6(p ~ 



2tt' 



One can perform the integration in Eq. (ITTI) and get 

{CuXkCk) ^ xPl'\P^ - P^) , (13) 

where x = S(ki/Cki denotes the fraction of extra fluctuation from the multi-stream effect. 

Now we compare the above result with the local form non-Gaussianity to estimate f^i 
in the multi-stream model. The local form f^i is defined as 

c = c. + ^/^l(cJ-(0), (14) 

or we can rewrite the above equation as 

((3) 

fNL ^ -(^2y^ ' (15) 
Comparing Eqs. (|T3|) and (|T5l) . we have 

fNL ^ XP-'" {^^^^^ ■ (16) 

Note that P^^^'^ ~ 10^ is a large number. We can have f^L — 100 when {P^ — P^)/P^ ~ 
3 X 10~^ and x ~ 0.3. As we shall show, this needs a tuning of the shape of the potential of 
order 1%. For f^L — 10, the tuning is of order 10%. Note that this non-Gaussianity appears 
only when one wave number is near ki ("near" denotes about one e-fold during inflation), 
fj^L is strongly scale dependent, and should have large running. 

We have not calculated the shape of non-Gaussianity for multi-stream inflation in the 
present paper. However, one note that non-Gaussianity is produced when k ^ ki, known as 
the squeezed limit in the literature. So we expect that the shape of non-Gaussianity should 
be similar to the local shape, with large running. So the local form estimator f^L should 
apply for our model. 

Finally, let us discuss the sign of f^i. In multi-stream inflation, the sign of /at^ is 
determined so that a positive 6(ki results in a larger perturbation at scale ki < k < k2- 
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\ . .fNL<0 

ax*;:: tiL>a 

* • . . B. . • • • ' ' ix 
► 

FIG. 2: We illustrate the scalar power spectrum as a function of scale k. The red (dashed) line 
denotes the trajectory A, and the blue (dotted) line denotes the trajectory B. If a perturbation 
6Cki > leads to path A, then we have f nl < on larger length scales and /nl > on smaller 
length scales. 

For example, when a fluctuation 6Cki > decides that the field roll to trajectory A, with 
Ca> Cb, then f^L > 0. 

It is also interesting to note that in the multi-stream model, J'nl can usually cross zero. 
It is because may, say, larger than on some scales, than become smaller than P^ on 
some other scales. As an example, when the deformation of the potential is "symmetric" , 

/ 6Vd^ = , (17) 

and the slow roll parameters satisfy 0{e) > 0{rf,^'^). In the lowest order of slow roll 
approximation, one have 

/ {P^ - Pf) = . (18) 

So in this example, Jnl must change sign at a certain scale. This special kind of running of 
the non-Gaussianity is illustrated in Figj2] 

IV. SUMMARY OF SIGNATURES FOR MULTI-STREAM INFLATION 

In this section, we will give a summary of signatures for multi-stream inflation. As we 
have mentioned in Section 2, 6(ki appears as an additional power at scale ki in the CMB. 
This additional power can not be arbitrarily large, because otherwise it should have been 
observed in the experiments. To be in accordance with CMB experiments, we require this 
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anomaly power to satisfy bC^k^ < XmP(^ , where Xm is determined by the error bar of CMB 
experiments, which varies when we choose ki to be different scales. Usually, one needs some 
fine tunning to satisfy the above constraint. In this section, we also discuss this tuning issue 
in more detail. 

First of all, the amount of tuning depends on the detail of the potential, thus depends 
on the detail of string landscape. For example, if the trajectories A and B are completely 
symmetric, then there is no tuning at all. However, in this case, there are no new signatures 
for multi-stream inflation either. 

To analyze the relation between tuning and signatures for multi-stream inflation further, 
let us calculate the relation between 6N and 6V. We use the slow roll parameters 

We first investigate the parameter regime where e is not too small. To be more explicit, 
0{e) > 0(?7^,.^^), and with "symmetric" deformation of potential: 

/ SVd^ = . (20) 

With these assumptions, 6N can be written as 

SM fV + dV V \ r,V fSV SV'\ 

The integration in Eq. (|2T|) can be expanded in powers of the slow roll parameters. So the 
natural value of 6N is 

SN^.tur.i ^ (iVi - N^f (^-4r^ + V " 7) ^ • ^^^^ 

We find that SN is suppressed by slow roll parameters. Note that we have assumed / SV{ip — 
(pi)d{p ~ {ip2 — ipiYSV. If for some reasons we have / 6V{lp — ipi)dip <^ {ip2 — (piY6V, then 
SN can be much smaller than given in Eq. (!22l) . 

To summarize the experimental signatures for multi-stream inflation, we propose two 
inequalities to divide the parameter space of multi-stream inflation into four regions. 

As discussed at the beginning of this section, current CMB experiments do not allow large 

1 /2 

features. A condition 6(ki < XmP^ is needed for multi-stream inflation to be consistent 
with CMB observations. If the above inequality is nearly saturated, then CMB experiments 
in the near future can detect the CMB features at scale ki. 
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Another inequality comes from the comparison of SV jV and P^-^"^ /x, where x = 6C,k-il P^^"^ ■ 
As shown in Eqs. (j3]) and (fT6|) . this comparison determines whether /tvl is larger or smaller 
than thus detectable or not in the near future. 

There is no estimator as sharp as f^L for the power asymmetry. So before explicitly 
fitting the multi-stream inflation model with CMB data, it is less clear how large power 
asymmetry can be detected in the next generation of experiments. However, one should 
note that the detection of f^L in Eq- (fT6|) indicates a difference between P^ and P^ . 
So here, we take the existence of large non-Gaussianity to be the condition for observable 
amount of power asymmetry. One should note that the detection of power asymmetry is in 
principle easier than the detection of /at^, because the detection of power asymmetry does 
not need information for Ck^ . We are glad to see whether improved methods can be proposed 
for detecting the power asymmetry. 

Using the above two inequalities, the four parameter regions for multi-stream inflation 
are: 

• When 5Ck, < x„,pI'^ and 5V/V < pI^'^/x, there is no observable signature for multi- 
stream inflation. There may be an isocurvature perturbation, which is different from 
single field inflation. However, this isocurvature perturbation can also arise in simpler 
two field models. 

1 /2 1 /2 

• When 8C,kx ^ XmPc^ and 5V/V < P^ /x, there is one signature for multi-stream 
inflation: a feature in the scalar power spectrum at scale ki. Such features in the 
scalar power spectrum may be responsible for the features in the WMAP data. 

1 /2 1 /2 

• When S(ki ^ XmP^ and 5V/V > Pq /x, there are three signatures for multi-stream 
inflation: The first signature is as that in the above case: the feature in the spectrum 
at scale ki. The second signature is that there will be an asymmetry in the powers 
between ki < k < k2 in one patch of the sky of scale ki and another. This could 
be responsible for the hemispherical power asymmetry. The third prediction is non- 
Gaussianity. The non-Gaussianity is of order 

f.,^xP-'/'^-^>l. (23) 

As discussed in the former section, J'nl in our multi-stream model is largest in the 
squeezed limit because ki <^ k, and /atl usually crosses zero. 
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Now we estimate whether this case is natural or not. The natural value of SN can be 
estimated using Eq. fl22|) . 

5iV,,turai ^ (iVi - N2fO{e)Pl/^fNL/x , (24) 

In order that perturbation on the scale ki is not ruled out by experiments, we need 

C 



6(ki < XmP}^"^ ■ The amount of tuning for this parameter region is measured by 



i5iV„.wr.i i5iV„...„.i (iVi - iV2)'0(«)/jVL 
For example, if we have multi-stream inflation for 3 e-folds, e = (9(0.01), and = 0.1, 
we have ^Cfci /(^^naturai — I/Za^l- To be more explicit, there is one potential in 10 
potentials that can produce fNL = 10, and one in 100 that can produce /nl = 100. 
Assume that future experiments show fNL = 10, then 3 features in the scalar power 
spectrum can naturally achieve that. 

1 /2 1 /2 

• When 6(ki XmP(^ and 5V/V > -P^ /a;, there is large non-Gaussianity, asymmetry 
on powers in different patches of the sky, and no feature in the power spectrum. 
However, more fine tuning is needed for this possibility than the former one. 

There is potentially one more signature of multi-stream inflation: the correlation between 
the curvature perturbation and the isocurvature perturbation. However, as we need to go 
to detailed models to calculate this correlation, we shall not consider it in the present work. 
We hope to address this issue in the future. 

If e is even smaller than yf and y^, then we are unable to perform the expansion (|22l) 
in the slow roll parameters. In this case, 6N/N ~ 6V/V, and we need more fine tuning 
to get the desired non-Gaussianity and anisotropy for perturbation powers for multi-stream 
inflation. 



V. CONCLUSION 



To conclude, we have described the multi-stream inflation model, and calculated the 
density perturbation and non-Gaussianity of the model. 

We find that the model can be parameterized by the 6V/V and 5N. We show that signa- 
tures such as non-Gaussianity, features in the CMB and the hemispherical power asymmetry 
can be produced in corresponding parameter regions. 
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We also note that isocurvature perturbation at ipi and may produce interesting physics. 
We hope we can address these issues in some future work. We also hope to build more explicit 
models of multi-stream inflation, and investigate the signatures of multi-stream inflation in 
more detail in the future. Finally, the perturbations are studied semi-quantitatively in this 
paper. It is important to calculate the perturbations more precisely and determine the 0{1) 
coefficients. To do so, one needs to investigate the cosmic perturbation theory for a quantum 
state with two decohered branches. We also leave this calculation to future work. 
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